Universal c-axis conductivity of high-T c oxides in the superconducting state 
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The anisotropy in the temperature dependence of the in-plane and c-axis conductivities of high-T c 
cuprates in the superconducting state is shown to be consistent with a strong in-plane momentum 
dependence of both the quasiparticle scattering rate and the interlayer hopping integral. Applying 
the cold spot scattering model recently proposed by Ioffe and Millis to the superconducting state, 
we find that the c-axis conductivity varies approximately as T 3 in an intermediate temperature 
regime, in good agreement with the experimental result for optimally doped YBa2Cu307-a; and 
Bi 2 Sr 2 CaCu208 + x. 



Microwave surface impedance measurements have pro- 
vided important information on the pairing symmetry 
|jj and quasiparticle relaxation [g| in the superconduct- 
ing state of high-T c oxides. For YBa 2 Cu 3 07_ a; (YBCO) 
(and similarly for other high-T c compounds the in- 
plane microwave conductivity a a b exhibits a large peak 
centered at approximately 25K in the superconducting 
state This peak structure of a a b is due to a com- 

petition between the quasiparticles life time and the nor- 
mal fluid density. From T c down to 25K the quasi- 
particle life time increases much more rapidly than the 
decrease of normal fluid density, causing a a b to rise 
with decreasing temperature. At low temperature, the 
quasiparticle life time reaches a limit and increases very 
slowly but the normal fluid density continues to fall, 
(Tab therefore falls with decreasing temperature. Along 
the c-axis, the conductivity behaves very differently ||. 
In contrast to its in-plane counterpart, a c falls below 
T c and does not show a conductivity peak. In opti- 
mally doped YBCO, a c rises slightly below 20K. But 
in Ba 2 Sr 2 CaCu 2 8 (BSCCO) 0,|^La 2 _ x Sr x CuO 4 [§, 
Tl 2 Ba 2 Cu0 6 and Tl 2 Ba 2 CaCu 2 8 (lOj, a c drops contin- 
uously from T c down to very low temperature, and no 
upturn appears. 

The different temperature dependences of a a b and a c 
are not what one might expect within the conventional 
theory of anisotropic superconductors. In this paper we 
present a detailed theoretical analysis for the c-axis con- 
ductivity. We shall show that the decrease of a c imme- 
diately below T c is explained by the fact that the region 
near the nodes, where the long lived quasiparticles exist, 
does not enter into the c-axis transport because of the 
anisotropic interlayer hopping integral. 

Let us first consider the behaviour of quasiparticle scat- 
tering in high-Tc materials. An important feature re- 
vealed by photoemission measurements is that the life 
time of quasiparticles is long along the Brillouin-zone di- 
agonals and short along other directions on the Fermi sur- 
face [O in both the normal and supercondcuting states. 



Based on this experimental result and the anisotropic 
temperature dependence of the in-plane and c-axis re- 
sistivity, Ioffe and Millis (IM) @,||] proposed a cold 
spot model to account for the normal state transport 
data. They assumed that the scattering rate of quasi- 
particles contains a large angular dependent part that 
vanishes quadratically as the momentum approaches the 
(0, 0) — (tt, 7t) line with negligible frequency and tempera- 
ture dependence and an isotropic but temperature depen- 
dent part, i.e. Tg = T cos 2 26> + r _1 , where 9 is the angle 
between the in-plane momentum of the electron and the 
a-axis. This type of the scattering rate was also used by 
Hussey et al. in the analysis of the angular dependent c- 
axis magnetoresistance . In Rcf . |Q , Ioffe and Millis 
has further assumed that r _1 has the conventional Fermi 
liquid form Tp L = j^+T~^ p . With this phenomenological 
model, they gave a good explanation for the temperature 
dependences of several transport coefficients in the nor- 
mal state. Van der Marel pl| has recently shown that 
this model also provides a good description for both the 
in-plane and c-axis optical conductivities in the normal 
state. 

The cold spot scattering rate, as discussed by IM jl2| , 
may be caused by interaction of electrons with nearly 
singular d x 2_ y 2 pairing fluctuations. In the supercon- 
ducting state, as the d x i_ y i -wave channel scattering is 
enhanced, the assumption of cold spot scattering made 
by IM is strengthened. This is indeed consistent with the 
recent photoemission data measured by Valla et al [ JLl| . 
Thus a detailed comparison between theoretical calcula- 
tions and experimental measurements for the transport 
coefficients in the superconducting state provides a cru- 
cial test for the cold spot model. 

In high-Tc superconductors, the electronic structure is 
highly anisotropic. In particular, the interlayer hopping 
integral depends strongly on the in plane momentum of 
electrons. For tetragonal compounds, the c-axis hopping 
integral is shown to have the form flq- |l9| , |2"2|| 
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t c = -t± cos 2 (20). 



(1) 



This anisotropic interlayer hopping integral is a basic 
property of high-T c materials. It results from the hy- 
bridization between the bonding O 2p orbitals and vir- 
tual Cu 4s orbitals in each CuC>2 plane and holds for 
all high-T c cuprates with tetragonal symmetry, inde- 
pendent of the number of CuC>2 layers per unit cell 
p8|^9[ . This form of the c-axis hopping integral was 
first found in the band structure calculations of high-T c 
oxides However, as shown in Refs. |l8|]l9| ], it is 

valid irrespective of the approximations used in these cal- 
culations. For Hg2BaCu04 or other non-body centered 
tetragonal compounds, t± is approximately independent 
of 9. For a body-centered tetragonal compound, such as 
La2- 2; Sr a ;Cu04, t± does depend on 9, but in the vicinity 
of the gap nodes it is finite. Since in the superconduct- 
ing state the physical properties are mainly determined 
by the quasiparticle excitations near the gap nodes, for 
simplicity we ignore the (9-dependence of t± in the dis- 
cussion given below. 

In YBCO, the CuO planes are dimpled with displace- 
ments of O in the c direction. O displacements, together 
with the CuO chains in YBCO, reduce the crystal sym- 
metry and introduce a finite hybridization between the a 
and 7r bands. This hybridization results in a small but fi- 
nite t c along zone diagonals which will change the low 
temperature behavior of the electromagnetic response 
functions. However, at not too low temperatures, Eq. 
(|l|) is still a good approximation. 

The conductivity is determined by the imaginary part 
of the current-current correlation function. If vertex cor- 
rections are ignored (23), it can be shown that the con- 
ductivity is given by 
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M(9) = ^Re A 2 cos 2 2^ 

' " \ia + iTgf - A 2 cos 2 29 1 
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where u ab {9) = 1, u c {9) = cos 2 (2(9), a ab = e 2 v 2 F N{0)/A, 
a c = e 2 t 2 ± N(0)/A, 7V(0) is the density of states of elec- 
trons at the Fermi level, and f(w) is the Fermi function. 
In obtaining the above equation, the retarded Green's 
function of the electron, G re t(k,uj), is assumed to be 



G ret (k,oj) 
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^ - ffcT 3 - AflTi + iTg 
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where = e a b{k) — t± cos k z u c {9) is the energy disper- 
sion of the electron, Ag — Ao cos 29 is the d-wave gap 
parameter, Tg is the quasiparticle scattering rate, T\ and 
r 3 are the Pauli matrices. 

In the normal state, A = 0, and Eq. (g) becomes 
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If the scattering is isotropic, i.e. Tg = T(T) independent 
of 9, then a a b and a c should have the same temperature 
dependence, in contradiction with experiments. How- 
ever, if Tg = T cos 2 29 + r -1 (T), then 



Cab = 



OLab 



2T- 1 2T" 1 




r n 



To 



(6) 
(7) 



When r t , a a b is proportional to yfr not r. This is 
the result first obtained by IM with a Boltzman equation 
analysis. If r _1 varies quadratically with T as in conven- 
tional Fermi liquid theory, the resistivity, i.e. a~ b , varies 
linearly with T. This provides a phenomenological ac- 
count for the linear resistivity of optimal doped cuprates. 
a c depends on two parameters, a c /To and tTq. In the 
limit Tq ^> r , a c depends very weakly on T and ex- 
trapolates to a finite value a c / (2Tq) at T = OK, in qual- 
itative agreement with the experimental data. These re- 
sults indicate that the simple cold spot model captures 
the key features of high-T" c transport properties in the 
normal state, although its microscopic mechanism is still 
unclear. 

In the superconducting state, Eq. (||) cannot be in- 
tegrated out analytically. However, in the temperature 
regime T c > T ^> Tq , where to = (r^ 1 ) is the thermal 
average of Tq, the leading order approximation in Tg is 
valid and the conductivity is given by [E3| 
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Tq 1 can be estimated from the experimental data of 
the in-plane microwave conductivity a a b and the nor- 
mal fluid density with the generalized Drude formula [|| : 
&ab ~ n a bT . For optimally doped YBCO ||, t^ 1 is less 
than IK at low temperatures and increases with increas- 
ing temperatures. At 60K, Tq" 1 is about 6K. Close to T c , 
Tq 1 becomes larger but still much less than the temper- 
ature. This means that the leading order approximation 
in Tg is valid in nearly the whole temperature range in 
which the experimental measurements have been done so 
far, at least for optimally doped YBCO. 

If Tg = T(T) = r-^T) does not depend on 9, Eq. (§) 
can be simplified to 



2 MT)t 
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where n^T) is the normal fluid density which decreases 
with decreasing temperatures. This is nothing but the 
generalized Drude formula which was first used by Bonn 
et al. in their data analysis for the in-planc microwave 
conductivity in the superconducting state |23|,|24f . From 



2 



Eq. it is easy to show that the ratio of the in- 

and out-of-plane conductivities <J a b/&c is proportional to 
the ratio n a b/n Cl i.e. cr a b/c c = n a b/n c . However, this 
does not agree with experiments, even qualitatively. It 
implies that the scattering rate must be anisotropic, as 
mentioned previously. 

In the cold spot model, T e = T cos 2 29 + t _1 (T) , 
(Jab and <t c behave very differently. Eq. (|J) now can be 
approximately written as 
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where £(3) = 1.202. In the high temperature limit 
T rT 2 /A 2 » 1, 



(12) 
(13) 



Gab ~ &ab\j -w-, j 
L 



where er„, c (01f) = a c /2r is the extrapolated normal 
state c-axis conductivity at Oif. These equations reveal 
a few interesting properties of the conductivities. Firstly, 
a ao is proportional to yfr and does not depend explicitly 
on Ao. This yfr dependence of a a b is an extension of 
Eq. @ in the superconducting state, which means that 
a ah (excluding the fluctuation peak at T c ) will change 
smoothly across T c since r changes continuously at T c . 
The temperature dependence of r in the superconduct- 
ing state is unknown. But phenomenologically it can be 
determined from the measured in-plane conductivity via 
Eq. flU). Secondly <7 C decreases monotonically with de- 
creasing temperature and behaves approximately as T 3 in 
the above temperature regime (Ao depends very weakly 
on temperature except close to T c ) . Furthermore a c does 
not depend on r, which means that this T 3 behavior is 
universal and independent of the impurity scattering pro- 
vided it is sufficiently weak that the coherent interlayer 
tunneling dominates. Furthermore, there is no free ad- 
justable parameters in Eq. ( O ) since both cr„ jC (0-ftf) and 
Ao can be determined directly from experiments. This 
therefore provides a good opportunity to test the cold 
spot scattering model by comparison with experiments. 

In the low temperature limit T tT 2 /AI < 1, Eq. (|) 
leads to the following results 
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In this limit both a a b and a c do not depend on To. This 
means that the To cos 2 29 term in Tg is not important in 
this temperature regime. In fact, Eqs. dlj) and (Il5|) are 



just the results of o a b and a c in an isotropic scattering 
system as given by Eq. (^|) since the normal fluid densi- 
ties n a b and n c behave as T and T 5 for tetragonal high-T c 
compounds at low temperatures, respectively fl^-|21|. In 
real materials where impurity scattering is not negligible, 
as discussed in Refs. [ |18|Jl9[] , this T 5 behaviour of the c- 
axis normal fluid density is fairly unstable and will be 
replaced by a T 2 law at low temperatures. In this case, 
the temperature dependence of a c will also be changed. 



The condition r rT 2 /Ag > 
v/TnT. Since TrT 1 



can 



be written as 
the condition 



T/A » i/Vr\ 

T/A > 1/VrVr holds if T/A > 1/VTo^ is satis- 
fied. If we assume T ~ 0.15ey , which is the value used 
by IM in their analysis of normal state transport coeffi- 



cients 1131, and r n 



6K as given by the experimental 



data for YBCO at 60K, then 1/yToTt) is estimated to be 
about 0.06. Thus the condition r rT 2 /A§ > 1 holds at 
least when T/A > 0.06 for YBCO. Since A - 2T C at 
low temperatures, therefore Eqs. (|l^) and ([t3]) are valid 
when T/T c > 0.12 for optimally doped YBCO. 
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FIG. 1. The c-axis conductivity as a function of (T/T c ) 3 
for YBa2Cu306.95 measured at 22GHz [5]. 

To compare the above results with experiments, we 
plot in Figure 1 the experimental data of a c at 22GHz 
as a function of (T/T c ) 3 for YBa 2 Cu 3 6 . 95 §. From 
30K up to T c , a c exhibits a T 3 behavior within exper- 
imental error, in agreement with Eq. (Jl3|). We have 
fitted the experimental data with other power laws of 
T/T c in the same temperature range, but found none 
of them fit the experimental data as well as the T 3 
power law. For this material, there is no experimental 
data on the temperature dependence of a c in the nor- 
mal state. But just above T c , a c sw 6.3 x 10 4 il _1 rn _1 
Since the normal state a c depends very weakly on 
T at low T for optimally doped YBCO @, we can 
therefore approximately take this value of a c as the ex- 
trapolated normal state c-axis conductivity at OK, i.e. 
a n , c (QK) « 6.3 x lO 4 ^- 1 ™- 
into Eq. (|l3|), we obtain a c r: 



i . Substituting this value 
6.8 x 10 5 (T/A ) 3 tt- 1 m- 1 . 
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By fitting this theoretical result with the corresponding 
experimental data in Figure 1, we find that Aq/T c w 2.6. 
This value of A /T c agrees with all other published data 
for optimally doped YBCO within experimental uncer- 
tainty. This agreement indicates that not only the lead- 
ing temperature dependence but also the absolute values 
of <t c predicted by Eq. ( |l3| ) agrees with experiments. 

For YBCO, t c is small but finite at the gap nodes. This 
has not been considered in the above discussion and if in- 
cluded, the temperature dependence of a c at low temper- 
atures will be changed. What effect leads to the upturn 
of a c at low T remains a puzzle to us. The model pre- 
sented here is inadequate to address this issue. Perhaps 
the interplay between the CuO chains and CuC-2 is play- 
ing an important role. A possible cause for the upturn 
of <j c is the proximity effect between CuO chains and 
Cu02 planes. However, a detailed analysis for this is too 
complicated to carry out at present and at this point we 
are unable to show explicitly whether this effect alone 
can lead to the observed upturn in a c . Another possi- 
bility is that this upturn is due to the onset of coherent 
tunnelling of the c-axis hopping whereas before this up- 
turn the c-axis hopping is incoherent. If this is the case, 
the assumption made in this paper would fail. However, 
since the dramatic increase in the quasiparticle lifetime 
P,|lEa] occurs at a temperature much higher t han the 
a c upturn temperature, we believe that this possibility is 
unlikely to be relevant. 
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FIG. 2. The c-axis conductivity as a function of (T/T c ) 3 for 
Bi 2 Sr 2 CaCu2 08+x (T c = 78K) [7]. The solid line is provided 
as a guide only. 

Figure 2 shows the c-axis quasiparticle conductivity 
obtained by Latyshev et al. from an intrinsic mesa tun- 
neling measurement for BSCCO. a c of BSCCO also 
exhibits a T 3 behavior in a broad temperature range 
in the superconducting state. However, the onset tem- 
perature of this T 3 term (~ 45K) is higher than that 
for YBCO. This is because BSCCO is very anisotropic 



and disorder effects are stronger than in YBCO. The 
crossover temperature from the impurity dominated limit 
at low temperatures to the intrinsic limit at high tem- 
peratures, as estimated in Ref. Jj, is about 30K. The 
value of a n>c (0K) for this material is difficult to deter- 
mine because of the pseudogap effect. Since the open- 
ing of a pseudogap always reduces the value of a ritC in 
the normal state, we can therefore take the value of 
&n,c just above T c (~ 3£l~ 1 m~ 1 ), as a lower bound to 
a ntC (0K). The c-axis conductivity measured at a voltage 
well above the pseudogap is nearly temperature indepen- 
dent and higher than the corresponding conductivity in 
the limit V -> (Figure 2 of Ref. 0). This conduc- 
tivity, o n ,c{eV > Ao) ~ 8f2 -1 m -1 , sets a upper bound 
to a n c (0K). Thus (T n , c (0K) is between 30 — 1 ttt,~ 1 and 
8ft m . By fitting the experimental data of a c from 
45K to T c with Eq. [t| we find that A /T c is within (2.3, 
3.1). This range of Ao/T c is consistent with the published 
data for BSCCO within experimental uncertainty. 

In conclusion, we have studied the temperature be- 
havior of microwave conductivity of high-T c cuprates in 
the superconducting state within the framework of low 
energy electromagnetic response theory of superconduct- 
ing quasiparticles. We found that the c-axis conductivity 
varies approximately as T 3 and does not depend on t(T) 
in the temperature regime T c ^> T ^> Ao/yTor. This 
universal tempearture dependence of er c agrees quan- 
titatively with the experimental data for YBCO and 
BSCCO. Our study shows that it is important to in- 
clude the anisotropy of the interlayer hopping integral 
in the analysis of c-axis transport properties of high-T c 
cuprates. 
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